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OF Z2 LATTICE MODELS IN THE NON-UNIQUENESS REGION 
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Abstract. We develop a new pressure representation theorem for nearest-neighbour 
Gibbs interactions and apply this to obtain the existence of efficient algorithms for ap¬ 
proximating the pressure in the 2-dimensional ferromagnetic Potts, multi-type Widom- 
Rowlinson and hard-core models. For Potts, our results apply to every inverse temperature 
but the critical. For Widom-Rowlinson and hard-core, they apply to certain subsets of both 
the subcritical and supercritical regions. The main novelty of our work is in the latter. 
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1. Introduction 


The pressure of an interaction is a crucial quantity studied in statistical mechanics and 
dynamical systems. In the former, it coincides with the specific Gibbs free energy of a 
statistical mechanical system (e.g. [16, Part III] and [35, Chapter 3-4]). In the latter, it 
is a generalization of topological entropy and has many applications in a wide variety of 
classes of dynamical systems, ranging from symbolic to smooth systems (e.g. [8, 23, 40]). 

In this paper, we continue the development in [15, 29, 9] of representing pressure with 
a simplified expression and using this to prove the existence of efficient algorithms for 
approximating pressure. 

We consider nearest-neighbour (n.n.) real-valued interactions <I> on , i.e. interactions 
defined only on configurations on single sites and pairs of adjacent sites. Since pressure 
is normally defined for stationary interactions, we assume that our interactions are sta¬ 
tionary here. Also, we allow the possibility of forbidden configurations S' on pairs of 
adjacent sites, and so the space of feasible configurations on ifi may be constrained. In 
the dynamical systems literature, the space of such feasible configurations is known as a 
nearest-neighbour shift of finite type (n.n. SFT), that here we denote Q.{S) (see Section 
3.1). 

A specification n for a n.n. interaction <t> is a uniquely determined collection of Borel 
probability measures given in an explicit form in terms of O, for configurations on finite 
subsets A of ifi and feasible configurations ^ on the boundary of A. A Gibbs measure jX for 
a n.n. interaction <I> is a Borel probability measure on Q.{S), whose conditional probability 
distributions on any such A agree with the specification for <I> for all boundary conditions 
^ of positive /r-measure. 

Gibbs measures exist for all n.n. interactions (and, indeed, for much more general 
interactions), but a given n.n. interaction may have more than one Gibbs measure. In many 
cases, including the ones of most interest to us here, there is a n.n. interaction <t> which 
gives rise to a parameterized family of interactions uniqueness of Gibbs 

measures holds for sufficiently small ^ (the so-called subcritical region) and uniqueness 
fails for sufficiently large ^ (the so-called supercritical region). 

Given a n.n. interaction <I> on a n.n. SFT D.{S), we can associate an energy to any 
feasible configuration on a finite subset A of Z'^. The partition function Zj of <I> on A 
corresponds to the sum over all feasible configurations on A of a function (namely, e^') of 
their corresponding energy, and the pressure P(<I>) is defined as the asymptotic exponential 
growth rate of the partition function on an increasing sequence of boxes B„ which 
exhausts ifi, as « —?> oo. Note that ?(•!>) implicitly depends on Q.(S). 

When d = 1, there is a closed-form expression for P('I>) in terms of the largest eigen¬ 
value of an adjacency matrix formed from <I> (see [28, p. 99]). In contrast, when d >2, 
there are very few n.n. interactions <I> for which P)*!*) is known exactly. 

There is much work in the literature on numerical approximations of P(<I>), both for 
somewhat general <I> and somewhat specific <I> (see [4, 14]). In our paper, we take a theo¬ 
retical computer science point of view (see [26]): an algorithm for computing a real number 
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r is said to be poly-time if for every S N, the algorithm outputs an approximation to 
r, which is guaranteed to be accurate within ^ and takes time at most polynomial in N to 
compute. In that case, we say that r is poly-time computable. 

One of our goals is to prove the existence of poly-time algorithms for P('I>) under certain 
assumptions on <I> and While one might expect such algorithms to exist for most O 

and Q.{S’) of practical interest, there exist Q.{S’) for which even P(0) (which corresponds 
to the topological entropy of when the n.n. interaction is <I> = 0) is not poly-time 

computable and some for which P(0) is not computable at any rate (see [22]). However, the 
closed-form expression when d = I mentioned above, always gives a poly-time algorithm 
in that case. 

We follow an approach initiated by Gamarnik and Katz [15], and further developed by 
two of the authors [29] of the present paper. The basic idea is motivated by the variational 
principle [23, Section 4.4], which asserts that P(<t>) is the supremum over all stationary 
Borel probability measures fi on D.{S') of the sum of two quantities: one quantity is the 
measure-theoretic entropy h{jj.) of ji and the other quantity is the integral, with respect to 
ji, of a simple explicit function Ao : Q.{S‘) M, determined by O. The entropy li{ji) can 
be expressed as the integral, also with respect to /i, of a function known as the information 
function 7^, i.e. /j(/i) = /The supremum is always achieved by a Gibbs measure ji 
for <I>, and so for such ji, we can write P('I>) = / (7^ 4-A<i>)(7/i. 

The idea of [15] was to represent P(<I>) as the integral of the same integrand, but with 
respect to a simpler measure v, i.e. P(<h) = / (7^ -\-Aci,)dv. This is what we call a pressure 
representation and requires some assumptions on p, v and Q(rf). 

A pressure representation becomes especially useful for approximating P(<I>) in the case 
that V is a periodic point measure, i.e. a measure which assigns equal weight to each dis¬ 
tinct translation of a given periodic configuration (this was the only case considered in 
[15]). Then / (7^ +A^)dv becomes a hnite sum. The terms in this sum corresponding 
to A<j, are easy to compute. In this way, the problem of approximating P(<I>) (and there¬ 
fore proving that P('l>) is poly-time computable) reduces to approximating 7^ on a single 
periodic conhguration and its translates. 

The pressure representation theorems in [15] and [29], as well as in our paper (see 
Theorem 6.3), work in all dimensions d. Among other conditions, these results require 
conditions on Q(ff) and a convergence condition for certain sequences of finite volume 
half-plane measures (different convergence conditions in the different results). In the case 
i7 = 2, if the convergence holds at exponential rate, then one obtains a poly-time algorithm 
for approximating P(4>) (see Theorem 9.1). For d > 2, one can deduce an algorithm for 
approximating P(<I>) with sub-exponential but not polynomial rate. 

In [ 15] and [29], the convergence condition is given in terms of the information function 
7^ of a stationary Gibbs measure p for the interaction. In our paper, the condition is given 
in terms of a closely related function which depends only on the specification n of 
the interaction (see Section 6.2), in contrast with [15] and [29]. This is natural, since the 
pressure depends only on the interaction and not on any particular Gibbs measure p. 

In [15], the convergence condition is strong spatial mixing of a Gibbs measure p for the 
n.n. interaction <I>. This condition is known to imply that there is a unique Gibbs measure 
for <I> and thus can be applied only in the uniqueness (subcritical) region of a given model. 
The convergence conditions in [29] are weaker but also apply primarily to this region. 
However, in our paper, since our convergence condition depends only on the interaction, 
one might expect that the pressure representation and approximation results can apply in 
the non-uniqueness region as well. Indeed, they do. As illustrations, we apply these results 
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to explicit subcritical and supercritical sub-regions of the 2-dimensional (ferromagnetic) 
Potts, (multi-type) Widom-Rowlinson and hard-core models. In particular, for the pressure 
approximation results for these models, we establish the required exponential convergence 
conditions. However, we believe that our results are applicable to a much broader class 
of models, in particular satisfying weaker conditions on (e.g. the topological strong 
spatial mixing property, introduced in [9]). We remark that the strong spatial mixing con¬ 
dition of [15] is a much stronger version of our condition, and so in this sense our results 
generalize some results of that paper (in particular, for the hard-core model on 1?). 

In the case of the 2-dimensional ferromagnetic Potts model, we obtain a pressure repre¬ 
sentation and efficient pressure approximation for all j3 f Pdq), where q is the number of 
colours, p is the inverse temperature and Pdq) = log(l -f ^/q) is the critical value which 
separates the uniqueness and non-uniqueness regions. Our proof in the non-uniqueness 
region generalizes a result from [11] for q = 2 (i.e. the Ising model) and we closely follow 
their proof, which relies heavily on a coupling with the bond random-cluster model and 
planar duality. For the uniqueness region, our result follows from [3]. (See Corollary 2, 
part 1.) 

For the Widom-Rowlinson and hard-core models, our results are not as complete as in 
the Potts case, since the subcritical and supercritical regions for these two models haven’t 
been completely determined, in contrast with the Potts model. We also expect our re¬ 
sults can be improved, because they only apply to proper subsets of the currently known 
uniqueness/non-uniqueness regions. 

For the Widom-Rowlinson model, in the supercritical region, we use a variation of 
the disagreement percolation technique introduced in [7], combined with the connection 
between the Widom-Rowlinson model and the site random-cluster model. In the subcritical 
region, we apply directly the results in [7]. (See Corollary 2, part 2.) 

For the hard-core model, in the supercritical region, we combine the coupling in [7] and 
a Peierls argument used by Dobrushin (see [13]). In the subcritical region, we use a recent 
result on strong spatial mixing for the hard-core model in Z^. (See Corollary 2, part 3.) 

For the Potts model, we also extend the pressure representation, by a continuity ar¬ 
gument, to give an expression for the pressure at criticality. It is of interest that there is 
an exact, explicit, but non-rigorous, formula for the pressure at criticality due to Baxter 
[5]. So, our rigorously obtained expression should agree with that formula, though we do 
not know how to prove this statement. It seems that Baxter’s explicit expression gives a 
poly-time approximation algorithm, but we cannot justify that our expression is poly-time 
computable. 

We remark that the finite volume half-plane measures mentioned above typically are 
constant on their bottom boundaries and thus are related to wetting models (see [34, 38]). 
Our proofs are related with such models where the interaction with the hard-wall is the 
same as the bulk interaction. 

The remainder of the paper is organized as follows. Since we have drawn heavily on 
many concepts from many different sources, for the convenience of the reader we have 
collected a good deal of relevant background material early in the paper. This can be found 
in Section 2, Section 3, Section 4, Section 5 and Section 7, with the notable exception of 
Lemma 5.5 in Section 5, there is very little new material in those sections. In Section 2 
and Section 3, we review the fundamentals on configuration spaces on id, Gibbs measures 
and pressure. In Section 4, we review the specific lattice spin systems models to which we 
apply our main results, and in Section 5 we review the bond and site random-cluster mod¬ 
els which are intimately connected with two of our models. Our pressure representation 
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theorem is contained in Section 6. We review spatial mixing and stochastic dominance in 
Section 7 and use these concepts in Section 8 to help establish exponential convergence 
results for our models. Finally, in Section 9, we combine our pressure representation theo¬ 
rem and our exponential convergence results in Section 8 to obtain pressure representations 
and poly-time algorithms for our models. 

2. Definitions and preliminaries 

2.1. Hypercubic lattice Z'^. Given £ N, we consider the d-dimensional hypercubic 
lattice Z"', which can be regarded as a countable graph with regular degree 2d, where 
V[W/) = is the set of sites and E{l/) = : x,y £ Z"^, ||x —y|| = l} is the set of 

bonds, with ||x|| = |x,| the 1-norm. We will mainly focus our attention on the case 

d = 2. 

Two sites x,y are adjacent if {x,y} £ E{!/) and we will denote this by x ~ y. All 
subsets of sites in will be denoted with uppercase Greek letters (e.g. A, A, 0, etc.). 
Whenever a hnite set A is contained in an inhnite set A, we denote this by A (g A. The 
(outer) boundary of A C Z^ is the set dA of x £ Z"' \ A which are adjacent to some element 
of A, i.e. 5A := {x £ A'^: dist({x}, A) = 1}, where dist(Ai,A 2 ) = niin;(gAj 
for Ai,A 2 C I/. We also write the closure o/A as A := AU (9A. On the other hand, 
the inner boundary of A C Z"^ is the set 5A (9A^ of x £ A which are adjacent to some 
element of A'^. When denoting subsets of 1/ that are singletons, brackets will be usually 
omitted, e.g. dist(x, A) will be regarded to be the same as dist({x}, A). 

A path T g Z'^ will be any sequence of distinct sites xi,... ,x„ such that x,- ~ x,+i, for 
all 1 < i < n. Similarly, a circuit C g Z'^ will be any path xi,... ,Xn with n > 4 such that, 
in addition, x„ ^ xi. We will say that the circuit is simple if x,- ^ Xj iff \i — J\ — 1 or 
{i,J} = {l,n} (in particular, xi,... ,x„ are all distinct). For A,© C Z'^, a path from Ato& 
is a path T whose first site is in A and whose last site is in 0. A set A C Z'^ is said to be 
connected if for every x,y £ A, there is a path T fromx to y contained in A (i.e. T C A). A 
set A g Z^ is said to be simply lattice-connected if A and A'^ are both connected. 

In Z'^ we can also dehne an alternative notion of adjacency and therefore, an alternative 
notion of boundary, inner boundary, closure, path, connectedness, etc., by replacing the 
1-norm || • || with the oo-norm || • ||oo, dehned as ||x||o<, = max;=i £/ |x,j, for x £ Z'^. When 
referring to these notions with respect to the oo-norm, we will always add a * superscript 
and talk about *-adj acency x ~y, ^-boundary d*A, inner ^-boundary 5*A, ^-closure A , -k- 
path, ^-connectedness, etc. Notice that two sites x and y are ★-adjacent if they are adjacent 
in a version of the c/-dimensional hypercubic lattice Z'^ including in addition diagonal 
bonds. We will denote this version of the lattice by l/’*. 

A natural order on Z'^ is the so-called lexicographic order, where y ^ x (or x g y) if and 
only if y 7 ^ X and, for the smallest i for which y, x,, y, is strictly smaller than x,. We also 
denote y ^ x (or x ^ y) if y ^ x or y = x. Considering this order, we dehne the family of 
sets g Z"^ as; 

(2.1) := {x )= 0 : —y < X < z}, 

where y,z £ Z'^ are such that y,z > 0 (here 0 denotes the vector (0,...,0) £ Z'^ and >, 
the coordinate-wise comparison of vectors). In addition, given n £ N, we dehne the n- 
block as the set B,, := [—n,nY D Z"' and we abbreviate by the set Si„ i„ = B,i \ 
where {x£Z'^:x^0} denotes the (lexicographic) past of 1/ and 1, the vector 

(1,...,1)£Z^. 
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2.2. Configuration spaces. Consider a finite set of symbols si called the alphabet. A 
configuration is a map 9 : A s/, for some (d ACZ‘^ (i.e. 9 G s/^), which will be 
usually denoted with lowercase Greek letters 9, T, V. The set A is called the shape of 9, 
and a configuration will be said to be finite if its shape is finite. For any configuration 9 
with shape A and A C A, 0(A) denotes the restriction of 0 to A, i.e. the sub-configuration 
of 0 occupying A. We will usually save the Greek letters ^ and Tj to denote configurations 
whose shape is the boundary dA of some given set A. For Ai and A 2 disjoint sets, 0 G 
si^^ and T G 0T will be the configuration on Ai U A 2 defined by (0 t)(Ai) = 0 and 
( 0 t)(A 2 ) = T. For a G si and A C a^ denotes the configuration of all a’s on A. A 
point is a configuration with shape Z'^, i.e. an element of s/'^ , usually denoted with the 
Greek letter co. 

Given sets Ai,A 2 C Z‘', A C Ai nA 2 and a pair of configurations 0 G si^^, T G si^'^, 
we define the set of A-disagreement as; 

( 2 . 2 ) EA( 0 ,T):={.reA: 0 (a)^t(a)}, 

i.e. the set of sites in A where 0 and 77 differ. 

The map O’ : Z'^ x — >■ will be the shift action on s/'^'^ defined by (x, (o) 

(Jxico), where x G Z'^ and co G s/^‘’, with (Oj:(a))) (y) = a)(x + y), for y G Z‘^. We also 
extend the shift action o^ to configurations with arbitrary shapes, i.e. given 0 G si^, we 
define 0^(0) G as the configuration such that (Oj:(0)) (y) = 0(x+y), fory G A —x. 

Given a point o) G si'^ , we define its orbit as the set O{co) := {(yx{(o)}xeZ‘‘- 
say that a point co is periodic if | 0 (a))| < 


2.3. Borel probability measures. Given a configuration 0 G si^, we define the cylinder 
set [0]a ;= {co G : ©(A) = 0} (or just [0], if A is understood). We denote by the 
O-algebra generated by all the cylinder sets with shape A and set := ^id. 

A Borel probability measure /X on is a measure determined by its values on cylinder 
sets of finite configurations such that jj.{s/^ ) = 1 - Given a cylinder set [0], we will just 
write /x(0) for the value of /x([0]). The support of such a measure ji is defined as; 

(2.3) supp(^) ;= I® G si^‘‘ : li{co{A)) > 0, for all A g Z^'j . 

Given A CACZ‘‘ and a measure p on we denote by /x|^ the restriction (or pro¬ 
jection or marginalization) of p to iPts- 

A measure p is shift-invariant (or stationary) if p(ax(A)) = p(A), for all measurable 
sets A G and xGZ'^ . Given any point co G s/^‘‘ and A G , we define the delta-measure 
supported on CO as the measure; 


(2.4) 


5(u(A) 


1 if©€A, 

0 otherwise. 


If © is a periodic point with orbit 0(©) = {© 1 ,..., ©*.}, we define v® to be the shift- 
invariant Borel probability measure supported on 0(©) given by; 


(2.5) 


2.4. Markov random fields. 
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Definition 2.1. Given A C Z^, a probability measure p on is a Markov random field 
(A-MRF) if, for any subset © (g A, any 9 € any AgAs.f. (90nACACA\©, and 
any T £ with p (t) > 0 , it is the case that: 

(2.6) p(0|t) = p(0|T(5©nA)). 

In other words, an MRF is a measure where every finite configuration conditioned to its 
boundary is independent of the configuration on the complement. 

3. Specifications, Gibbs measures and pressure 

3.1. Gibbs specifications. Fix a dimension G N and let S = {S\,... ,Sd) be a set of 
constraints such that C for i = 1,... ,d. Given any set A C Z'^ and a configuration 
0 G.!^/^ we say that 0 is feasible for S’ if for every x G A such that {x,x + e, } C A, we have 
that (0(x), 0(x + e,)) ^ Si, where ei,... , 6 ^/ is the canonical basis. The nearest-neighbour 
shift of finite type (n.n. SFT) Q(S) induced by S, is the set of points: 

(3.1) i2(S) := |a) G : co is feasiblej . 

We will always assume that Q(S) f- 0. 

In the symbolic dynamics literature, a feasible configuration on a set A is called locally 
admissible, and is called globally admissible if it also extends to a point of i2(S). 

Notice that Q(S) is always a shift-invariant set, i.e. (Tx(i2(S)) = i2(S), for all x G 
Given a n.n. SFT i2(S), .M\ (i2(S)) denotes the set of Borel probability measures whose 
support supp(p) is contained in Q(S) and o-(n(<o’)) C .y£'i(Q(S)), the correspond¬ 
ing subset of shift-invariant Borel probability measures. Given a configuration 0 G 
j 0 ]n((S’) j.jjg ggj. (or just [0]^(‘^) if A is understood). 

Definition 3.1. A nearest-neighbour (n.n.) interaction/or a set of constraints S is a real¬ 
valued shift-invariant function <I> from the set of configurations on sites x and feasible 
configurations on bonds {x,x-|-e,} to ]R,/orx G and i = 1,... ,d. Here, shift-invariance 
means that <I>((j;t(0)) = ^(0) for configurations 0 on sites and bonds, and for all x G 

Often in the literature a n.n. interaction is not required to be shift-invariant. Our as¬ 
sumption of shift-invariance on a n.n. interaction fits naturally with the shift-invariance of 
a n.n. SFT. Clearly, a n.n. interaction is defined by only finitely many numbers, namely 
the values of the interaction on configurations on { 0 } and bonds { 0 ,e,}, for i = l,...,d. 

We can view an interaction <1> as implicitly determining the constraints S, and hence 
Q(S), by the absence of S from the domain of <I>. Some authors incorporate the constraints 
by allowing the interaction to take the value -|-oo. 

Definition 3.2. Given a n.n. interaction ^ for a set of constraints S and a set A (£ Z'^, we 
define the energy function EJ : R as: 

(3.2) Ej(0):=^<I>(0(x))-ff L <I>(0({x,x + p,})), 

x€A 1=1 {A:,x+e,}CA 

where 0 is any feasible configuration in We define the partition function of A as: 

(3.3) Zj:= ^ exp(-Ej(0)), 
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(3.4) 


n^/\e):= 


and the following boundary-free probability measure on 

' ^ exp (—E^(0)) if 0 is feasible, 

0 otherwise. 

Analogously, for an arbitrary (O G we can take ^ = a)(5A) and consider: 

(3.5) 


■^A • — 


L exp(-E^(0^)), 

6\ 6^ feasible 


and then define the ^-boundary probability measure on 


(3.6) 


4(0): = 


-exp 


(-E^( 0 ^ 


if 0^ is feasible, 

0 otherwise. 

The collection Tt = {4}a,§ called a 7/ Gibbs specification/or the n.n. interaction 
<I>. For A C A and T G we marginalize as follows: 

r^l 

ees/^:8{A)=t 


(3.7) 


4(^) = L 4(®)- 


Notice that each 4 is an MRE on In addition, a Gibbs specification n as defined 

above is always stationary, in the sense that n^^^\ax{A)) = n^{A), for every A C We 
will usually think of the set of restrictions ^ implicit when considering a n.n. interaction 
•t*. Given a point co € D.{S'), we will abbreviate: 

(3.8) 

3.2. Gibbs measures. 


Definition 3.3. A nearest-neighbour (n.n.) Gibbs measure/or a n.n. interaction <I> is a 
measure p G (£2(if)) such that for any A (s Z'^ and CO G with ll{co{dA)) > 0, we 
have that > q and: 

(3.9) M(0|=^Ac)(®) = <(0)M-fl-s., 

for 0 G where {4}a,i^ stationary 1/ Gibbs specification for <I>. 

While our interactions and specifications are assumed to be shift-invariant, a Gibbs 
measure for such an interaction may or may not be stationary. The definition of n.n. Gibbs 
measure, shows that such a measure is an MRE. The definition is stated only for cylinder 
events [0] in A, but this is equivalent to the usual definition with general events A G 
instead. 

Every n.n. interaction <t> has at least one (stationary) n.n. Gibbs measure (special case 
of a general result in [35, Theorem 3.7 and Theorem 4.2]). Eor a single <I>, multiple Gibbs 
measures can exist. This phenomenon is usually called a phase transition. 


3.3. Pressure. Now we proceed to define the pressure of a n.n. interaction <I>. 

Definition 3.4. Given a n.n. interaction ^ for a set of restrictions S', the pressure of <I> is 
defined as: 
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Given « S N, we can also define an analogous version Zg of the partition function Zg , 
but over globally admissible configurations: 

(3.11) Zl:= £ exp(-E^^(0)). 

Notice that Zg^ < Zg^. The following result states that in the normalized limit, both 
quantities coincide. 

Theorem 3.1 ([35, Theorem 3.4], see also [14, Theorem 2.5]). Given a n.n. interaction <I> 
for a set of restrictions S’: 

(3.12) = 

The pressure is the main quantity of interest in this paper. Our goals are to find simple 
representations of pressure in terms of very special configurations and use this to develop 
efficient (in principle) algorithms to approximate the pressure. 


4. Main models: Potts, Widom-Rowlinson and hard-core 


In this section we introduce the three main families of lattice models studied in this 
paper. The first one will be the Potts model, which can be regarded as a generalization 
of the Ising model by considering more than two types of particles. The second one, the 
Widom-Rowlinson model, is also a multi-type particle system but with hard-core exclusion 
between particles of different type. The third one is the classical hard-core model. 


4.1. The (ferromagnetic) Potts model. Given d,q£N and j3 > 0, the (ferromagnetic) 
Potts model with q types and inverse temperature j3 is defined over the alphabet = 
{ 1 ,..., and given by the n.n. interaction: 


(4.1) 


0 ^( 0 ) 


—]3 if 0 (x) = 0 (x-|-e,), 

0 if 0 (x) 7 ^ 0 (x-|-e/), 


for 0 G X G Z'^, i = l,...,d, where the constraints Si are empty. The speci- 

fication a f induced by <Pp defines the (ferromagnetic) Potts model, where 

neighbouring sites preferably align to each other with the same type or “colour” from the 
alphabet 

A measure p G is called a Potts Gibbs measure for q types and inverse tem¬ 

perature j3 > 0 if it is a n.n. Gibbs measure for the specification above. 


Theorem 4.1 ([6]). For the J? (ferromagnetic) Potts model with q types and inverse tem¬ 
perature j3, there exists a critical inverse temperature Pdq) ■= log(l + y/q) such that 
uniqueness of Gibbs measures holds for j3 < ^c{q) and for p > pc{q) there is a phase 
transition. 


4.2. The (multi-type) Widom-Rowlinson model. Given G N and Z > 0, the Z^ 

Widom-Rowlinson model with q types and activity X is defined over the alphabet SSq = 
{0,1,... and given by the set of constraints S =(S\^... ^Sfi), where ^' = (0 G \ {0})^ : 
0(1) d (^(2)}. for all / = 1,... ,^f, and by the n.n. interaction for S over configurations on 
sites: 


(4.2) 




-log(Z) if 0 G {1,...,^}, 
0 if 0 = 0 , 
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where 6 € and x € Z'^. The specification a}k^ induced by <I>;t defines 

the (multi-type) Widom-Rowlinson model, where neighbouring sites are forced to align to 
each other with the same type or “colour” from the alphabet or with 0 . 

A measure ji € ) is called a Widom-Rowlinson Gibbs measure for q types and 

activity X > 0 if it is a n.n. Gibbs measure for the specification above. 

Theorem 4.2 ([36], see also [19]). For the 7? Widom-Rowlinson model with q types and 
activity X, uniqueness of Gibbs measures holds for sujficiently small X and there is a phase 
transition for sujficiently large X. 

4.3. The hard-core lattice gas model. Given 7 > 0, the Z'^ hard-core model with activity 
7 is defined over the alphabet {0,1}, and given by the set of constraints S’, where (ai = 
{( 1 , 1 )}, for all / = 1 ,.. .,d, and the the n.n. interaction for S over configurations on sites: 


(4.3) 


<Py{e) 


-log( 7 ) if 0 = l, 
0 if 0 = 0 , 


for 0 € {0,1}^’^^, X G Z'^. The specification induced by <t>y defines the 

hard-core model, where neighbouring sites cannot be both 1 . 

A measure ji G ({0,1}* ) is called a hard-core Gibbs measure for activity 7 > 0 if 
it is a n.n. Gibbs measure for the specification above. 


Theorem 4.3 ([17, Theorem 3.3]). For the 1? hard-core model with activity 7 , unique¬ 
ness of Gibbs measures holds for sujficiently small 7 and there is a phase transition for 
sujficiently large 7 . 


For both the Potts and Widom-Rowlinson models we will also distinguish a particular 
type of particle or colour in the alphabet. W.l.o.g., we can take the type q'ms^qOX SSq \ (Oj, 
respectively. Given this colour, we will denote by (Oq the fixed point q^‘‘. For the hard-core 
model, we will consider the two special points and given by: 


(4.4) 


and = ( 7 e| 



if Y^iXi is even, 
if is odd. 


5. Random-cluster models 

The Potts and Widom-Rowlinson models have interpretations in terms of a random- 
cluster representation. The Potts model is related to a random-cluster model on bonds (via 
the so-called Edwards-Sokal coupling), while the Widom-Rowlinson is naturally related to 
a random-cluster model on sites. 

Definition 5.1. A coupling of two probability measures pi on a finite set X and p 2 on a 
finite set Y, is a probability measure P on the setXxY such that, for any A<GX and B G_Y, 
we have that: 


(5.1) 


P(A X 7) =pi(A) andf‘{XxB) = p 2 {B). 
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5.1. The bond random-cluster model and the Potts model. We will make use of the 

bond random-cluster model. One of our main results. Part I of Theorem 8.1, is proven 
using arguments based on this model. This model is a two parameter family of dependent 
bond percolation models on a finite graph. We are mainly interested in finite subgraphs of 
I? and we describe the model with boundary conditions indexed by i = 0,1. 

Fix a finite simply lattice-connected set of sites A. Let E^{A) denote the set of bonds 
with both endpoints in A and E * (A) the set of bonds with at least one endpoint in A. We 
speak of a bond e as being open if w{e) = 1, and as being closed if w(e) = 0. 

Definition 5.2. Given a finite simply lattice-connected set A, and parameters p S [0,1] 
and q > 0, we define the free (i = Q} and wired (i = Ij bond random-cluster distributions 
on E'{A) (i = 0,1) as the measures each w S {0, assigns probability 

proportional to: 

(5.2) n 

[eeEi{A) J 

where #i(w) is the number of open bonds in w and k'^{w) and k\^{w) are the number of 
connected components (including isolated sites) in the graphs (A, {e G E^{A) : w{e) = 1}) 
and [1?,E^{1?\A)\J {e G E^ {A) : w(e) = 1}), respectively. 

Notice that when ^ = 1, we recover the ordinary Bernoulli bond percolation measure 
tj^p.A, while other choices of q lead to dependence between bonds. For given p and q, 
one can also define bond random-cluster measures on 1? as a limit of finite volume 
measures 0^ ^ ^ (i = 0,1). 

Theorem 5.1 ([17, Lemma 6.8]). For p G [0,1] and q GN, the limiting measures: 

(5.3) 0«=lim0«_^^, iG{0,l}, 

exist and are translation invariant, where {A„}„ is any increasing sequence of finite simply 
lattice-connected sets that exhausts 7?. 

General bond random-cluster measures on 7? can be defined using an analogue of the 
DLR condition [20, Definition 4.29]. For ^ > 1, there is a value pc{q) that delimits exactly 
the transition for existence of an infinite open cluster for these measures. It is known [20, 
p. 107] that for q> 1 and p < Pc{q), there is a unique such measure which we denote by 
0p ^ (characterized by the nonexistence of infinite open clusters), and that coincides with 

(jtjyq and 0p*,j in this region. It was recently proven (see [6]) that pc{q) = t+^’ s^ery 

q>i- 

Let p = I — e^^. The/ree Edwards-Sokal coupling (see [20]) is a coupling be¬ 

tween the boundary-free Potts measure and 0^**^ The wired Edwards-Sokal coupling 

A a coupling between ;rjA ^^d 0pp a- Notice that pc{q) = 1 — 

These couplings are measures on pairs of site configurations and corresponding bond 
configurations. The projection to site configurations is the boundary-free/Op-boundary 
Potts measure, and the projection to bond configurations is the free/wired bond random- 
cluster measure, respectively. 

Theorem 5.2 ([20, Theorem 1.13]). Let A be a finite simply lattice-connected set, q gN, 
and let p G [0,1] and > Obe such that p = I — e^^. Then: 
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(1) For w € {0, the conditional measure a ^ {'^}) 

tained by putting random colours on entire clusters ofw not connected with 7?\A 
(of which there are (w) — 1) and colour q on the clusters connected with 1? \ A. 
These colours are constant on given clusters, are independent between clusters, 
and the random ones are uniformly distributed on the set 

(2) For 9 £ the conditional measure a {0} x {0, on {0, 

is obtained as follows. Consider the extended configuration 9 = 9q^^ and an 
arbitrary bond e = {x,y} £ ^'(A). If9[x)f^ 0(y), we setw{e) =0. If9[x) = 9{y), 
we set: 


(5.4) 


^(e) = 


1 with probability p, 


0 otherwise, 

the values of different w{e) being (conditionally) independent random variables. 


The couplings can be used to relate probabilities and expectations for the Potts model 
to corresponding events and expectations in the associated bond random-cluster model. A 
main example is a relation between the two-point correlation function in the Potts model 
and the connectivity function in the bond random-cluster model [20, Theorem 1.16]. 

By considering a displaced version of I?, namely j\+7? (the dual lattice), we can 
define a notion of duality for bond configurations w. Notice that every bond e € E (Z^) (if 
we think of bonds as unitary vertical and horizontal straight segments) is intersected per¬ 
pendicularly by one and only one dual bond e* £ £(^ 1 -f Z^), so there is a clear correspon¬ 
dence between E{Z^) and E{^ \ +1?). We are mainly interested in wired bond random- 
cluster distributions on the set of sites B„ := [—n -f 1, «]^ fl Z^. Given « £ N, if we consider 
the set of bonds £* (B„), it is easy to check that there is a correspondence e^ e* between 
this set and the set of bonds from -fZ^ with both endpoints in [—n,n]^n (^1 +Z^), 
which can be identified with the set £’*'(B„). Then, given a bond configuration w £ £* (B„) 
we can associate a dual bond configuration w* £ £’'*(B„) such that w*{e*) = 0 if and only 
if w{e) = 1. 

Considering this, we have the corresponding equality: 


Proposition 5.3 ([20, Equation 6.12 and Theorem 6.13]). Given n G N, p G [0,1] and 
q G N.- 


(5.5) 


CoA'i 


= <’", 3 . (“■*)' 


for any bond configuration w G {0,1}^'(®"^, wherein = [—n-f l,n]^nZ^ and p* G [0,1] 
is the dual value of p, which is given by: 


(5.6) 


p* g(i-p) 

1-p* p 


The previous duality result can be generalized to more arbitrary shapes and it has also a 
counterpart from free-to-wired boundary conditions, instead of from wired-to-free. 

The unique fixed point of the map p^ p* defined by (5.6) is and, as mentioned 
above, is known to coincide with the critical point pc{q) for the existence of an infinite 
open cluster for the bond random-cluster model (see [6, Theorem 1]). It is easy to see that 
p>Pc{q) iff p* <pc{q). 
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5.2. The site random-cluster model and the Widom-Rowlinson model. In a similar 
fashion to the bond random-cluster model, we can perturb Bernoulli site percolation, where 
the probability measure is changed in favour of configurations with many (for q> 1) or few 
(for q <l) connected components. The resulting model is called the site random-cluster 
model. 

Definition 5.3. Given A d 7?, and parameters p G [0,1] and q >0, the wired site random- 
cluster measure probability measure on {0,1}^ which to each 0 G {0,1}''^ 

assigns probability proportional to: 

(5-7) = 

where A = #i(0) is the number of 1 ’i in 0 and K/^(9) is the number of connected 

components in {x G A.: d{x) = 1} that do not intersect 5 A. 

The/ree site random-cluster measure defined as in (5.7) by replacing ffA(0) by 

the total number of connected components in A. However, we will not require that measure 
in this work. In any case, taking q = 1 gives the ordinary Bernoulli site percolation a, 
while other choices of q lead to dependence between sites, similarly to the bond random- 
cluster model. 


Proposition 5.4. Given a set A (s 7? and parameters A > 0 and q G N, consider the 
Widom-Rowlinson with q types distribution and monochromatic boundary condition 
Now, let f : > {0,1}'^ be defined site-wise as: 


(5.8) 


(/(0))W 


0 if0{x)=0, 

1 if eix)fo, 


for 0 G and x G A, and let p = -j^. Then, = V^^^A’ ^here f*'^i\{‘) ■= 

{')) denotes the push-forward measure on {0,1}^. 


The requirement that Ka( ) does not count connected components that intersect the inner 
boundary of A in the site random-cluster model, corresponds to the fact that non 0 sites 
adjacent to the monochromatic boundary 0)^(5A) in the Widom-Rowlinson model must 
have the same colour q. 

For q = 2. Proposition 5.4 is proven in [21, Lemma 5.1 (ii)], and the proof extends easily 
for general q. Proposition 5.4 can be regarded as a coupling between and V^^^A’ 
because a push-forward measure can be naturally coupled with the original measure. 

It is important to notice that ^ is itself not an MRF: given sites on a simple circuit C, 
the inside and outside of C are generally not conditionally independent, because knowledge 
of sites outside C could cause connected components of 1 ’s in C to “amalgamate” into a 
single component, which would affect the conditional distribution of configurations inside 
C. The following lemma shows that in certain situations, when conditioning on a circuit C 
labeled entirely by 1 ’s, this kind of amalgamation does not occur. 


Lemma 5.5. Let % f @ G_ A <^7? be such that A'^ is connected. Take A 5*0 DA. 

Consider an event A G and a configuration T G {0,1}^, where E C A \ 0 . Then: 




(5.9) 
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Proof. W.l.o.g., we may assume that A is a cylinder event [0] with 0 G {0,1}® (by linear¬ 
ity) and H = A \ 0 (by taking weighted averages). 

Now, H = A \ 0 can be written as a disjoint union of *-connected components E = 
K\U ■ ■ - UKn. For every i, d*Ki C A'^ U ^ (in fact, d*Ki C A'^ U A). Since A'^ U ^ is 
connected and A is finite, for every site in d*Ki there is a path to infinity that does not 
intersect Ki. 

Then, by application of a result of Kesten (see [24, Lemma 2.23]), d*Ki is connected, 
for every i. In addition, we have that A = 0 U A U E and d*Ki C A^ U A. 

We claim that: 

(5.10) ka(u) = K-A(t)(0)Uo^) + ^K-A:,(r)(/:,)) = ka('u(0)Uo^) + k-i(t), 

i=\ 

for any u G {0,1}"^ such that t)(A) = U and v(L) = r. 



Figure 1. A 7t-connected 0 (in black), the set A = 5*0 n A (in dark 
grey) and A"^ (in light grey) for A = Sy^z- 


To see this, given such i), we exhibit a bijection r between the connected components 
of V that do not intersect 0A and the union of: (a) the connected components of V (0) UO^ 
that do not intersect 0A, and (b) the connected components of v(Ki) that do not intersect 
for all /; namely, if C C A is a connected component of v, then r is defined as follows: 


(5.11) 


r(C) 


cn^ ifcn^T^e, 

C ifCCE. 


In order to see that r is well-defined, note that if C intersects 0 and E, the set C fl 0 
is still connected thanks to the fact that d*Ki is connected and u(A) = 1'^. To see that r 
is onto, observe that if C' is a connected component of u(0)UO^, then there is a unique 
component C of u such that Cfl 0 = C', due again to the fact that d*Ki is connected. And 
r is clearly injective because two distinct connected components cannot intersect. 

Finally, we conclude from (5.10) that: 


(5-12) V/ij,A(0 I Ut) = ----- 
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(5.13) 

(5.14) 

as we wanted. 




;(^#,{0iA)^k-a{0iAoA\0*) 




□ 


Remark 1. We claim that ;/ 0 ^ 0 C A (s are such that A'^ is connected, 0 is -k- 
connected and 0 fl 5A ^ 0, then A'^ U 0 is connected, which is the main hypothesis of 
Lemma 5.5. This follows from the easy fact that the k-closiire of a k-connected set is 
connected. 


6. Pressure representation 

6.1. Variational principle. The variational principle states that the pressure of an inter¬ 
action has a variational characterization in terms of shift-invariant measures. We state the 
variational principle below for the case of an n.n. interaction <I> for a set of restrictions S’. 

Theorem 6.1 (Variational principle [23, 33, 35]). Given a n.n. interaction <i>for a set of 
restrictions S, we have that: 

(6.1) P(4>) = sup (h{p)+[Aciydfi), 
where: 

• A<j,(tu) := —<I>(a)(0)) —LSLi^(®({0:e,})),/or O) € Li{S), and 

• h{p) := lim„^oo /r(0)log(^(0)) is the measure-theoretic entropy of 

p, where OlogO = 0. 

In this case, the supremum is also always achieved (see [23, Section 4.2]) and any 
measure which achieves the supremum is called an equilibrium state for A<i>. So, if p is an 
equilibrium state, then; 

(6.2) P{<P)=h{p) +J A^dp. 

For a shift-invariant measure p and A g Z"^ \ {0}, define: 

(6.3) p^,A(t0) :=^(a)(0)|a)(A)), 

and let Pfi{(0) := lim„^oop^_B„n.^(®)> which exists p-a.s. [23, Theorem 3.1.10] by Levy’s 
zero-one law. In addition, let; 

(6.4) /^(®) :=-logp^(®), 

which is also defined p-a.s. and is usually called the information function. It is well-known 
(see [16, p. 318, Equation 15.18] or [27, Theorem 2.4, p. 283]) that for any shift-invariant 
measure p, h{p) = / I^dp. Therefore, if p is an equilibrium state for <I>, we can rewrite 
the preceding formula for P(<I>) as: 

(6.5) Pi^) = J {If,+A^)dp. 

So, the pressure can be represented as the integral of a function, determined by an 
equilibrium state p and 4>, with respect to p. 

In this section, we show that the pressure can be represented as the integral of a function 
similar to If, -l-A<j,, with respect to any invariant measure v, assuming some conditions. 
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This is useful for approximation of pressure when V is an atomic measure supported on a 
periodic configuration (see Section 9). 

One of the conditions involves the SFT Q.{S'). 

Definition 6.1. A n.n. SFT Q.{S‘) for a set of constraints S satisfies the square block D- 
condition if there exists a sequence of integers {r„}„>i such that ^ 0 as « —> and, for 
any finite set A g ® G T € 

( 6 . 6 ) ^ ^ ^ 0 . 

This condition is a strengthened version of the classical D-condition (see [35, Section 
4.1]) which guarantees that the set of Gibbs measures for <I> coincides with the set of 
equilibrium states for A<i). 

Definition 6.2. Given a set of restrictions S', the corresponding n.n. SFT £2(<o’) C 
and a G xS, we say that £1(S) has a safe symbol a if {a,b), {b,a) ^ Sj, for every b € s^, 
for all i = I,... ,d. 

It is easy to see that if Q.{S) has a safe symbol, then it satisfies the square block D- 
condition. For the sets of restrictions S in the Potts, Widom-Rowlinson and hard-core 
models, the corresponding n.n. SFT Q.{S) has a safe symbol in each case (any a € 

0 G SSq, and 0 G {0,1}, respectively), so Q(S) satisfies the square block D-condition for 
the three models. 

6.2. The function A and additional notation. Given a n.n. interaction 4> for a set of 

constraints S, we will define some useful functions from Q(S) to R. First, given 0 G A (g 
Z'^ and CO G Ti(S), we define: 

(6.7) ;rA(a)) := <(0(0) = a)(0)) = = co{0)). 

Recall that, for y,z G Z'^ such that y,z> 0, we have defined the set S'jj j as {x )= 0 : 
—y < X < z}. Now, given y,z > 0 and co G T1(S), define 7ry-(co) Tts^^ico) and, given 
n G N, abbreviate K„{co) := 7ri„.i,i(cu). Considering this, we also define the limit 7t{co) := 
\im„^oo7t„{(0), whenever it exists. If such limit exists, we will also denote I„{co) := 
-\og7t{(o). 

It is not difficult to prove that under some mixing assumptions over an MRF fi, namely 
the SSM property introduced in Definition 7.1 (see Section 7), and assuming that supp(/r) = 
T1(S), one has that the original information function coincides with in i2(S). In this 
sense, our definition provides a generalization of previous results (see [15]), where may 
not be even well-defined. 

Now, suppose we have a shift-invariant measure v such that supp(v) C i2(S). We say 
that: 

(6.8) lim 7ry,z(ftt) = A(co) uniformly over a> G supp(v), 

if for all e > 0, exists k G N such that for all Vy,z > Ik: 

(6.9) \7ty^^{co) — n{co)\ < e, forall to G supp(v). 

In addition, we introduce the following bound: 

(6.10) C;[(v) := inf{7rA(tu): 0 G A (g Z^,a) G supp(v)}. 

Lemma 6.2. Let n be a n.n. interaction ^ for a set of restrictions S, with % and i2(S) the 
corresponding specification and n.n. SFT. Then, ifLl{S) has a a safe symbol, we have that 
C;r(v) > 0, for any shift-invariant measure V such that supp(v) C Q,(S). 
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Proof. The proof is analogous and a particular case of [29, Proposition 2.17]. In that 
reference, under these assumptions, it is shown that c^(v) := inf{^(a)(0)|a)(A)) : A d 
1‘^\{0},COG supp(v)} > 0, for a given n.n. Gibbs measure ji for <I>. We leave it to the 
reader to verify that C;t(v) > (v), for any such fi. □ 

In fact, much weaker conditions than the existence of a safe symbol are sufficient for 
the result of Lemma 6.2 and also for having the square block D-condition. See, for ex¬ 
ample, the single-site fillability property [29] and the topological strong spatial mixing 
property [9]. Notice that, since the Potts, Widom-Rowlinson and hard-core models have a 
safe symbol, we have that C 7 i{v) > 0, for any shift-invariant v with supp(v) C 

6.3. Pressure representation theorem. Pressure representation results can be found in 
[29, Theorems 3.1 and 3.6]. Those results are not adequate for the application to the specific 
models we are considering in this paper. Instead we will use the following result, whose 
proof is adapted from the proof of [29, Theorem 3.1], as well as an idea of [29, Theorem 
3.6]. In contrast to the results of [29], our result makes assumptions on the specification 
rather than a Gibbs measure. 


Theorem 6.3. Let ^ be a n.n. interaction for a set of restrictions S and suppose that 
satisfies the square block D-condition. Let V be a shift-invariant measure such that 
supp(v) C n((o’) and Cj[{v) > 0. In addition, suppose that: 

(6.11) lim Ky^^co) = k[(o) uniformly over tu G supp(v). 

y,z-^’=° ■ ’ 

Then: 

(6.12) Pi^)= J {I^+A^)dv. 


Proof. Choose £ < 0 and L > 0 to be lower and upper bounds respectively on values of <I>. 
Given n G N, let r„ be as in the definition of the square block D-condition and consider the 
sets B„ and A„ := We begin by proving that: 

(6.13) ^(logZf,, -t-log<(®(B„)) +El{(o{B„))) ^ 0, 

uniformly in O) G H. For this, we will only use the square block D-condition. We fix n G N, 
(0 G supp(v) and let m„ := |A„| — |B„|. Let > 1 be a constant such that for any AdZ'', 
the total number of sites and bonds contained in A is bounded from above by Crf|A|. 


(6.14) <„(a)(B„))><,^(®(A„)) 

exp(-E^(®(A;))) 

15) __—-_ 

^G-.6co{dA„) feasible ®Xp( £—(0fi)((9A„))) 

j ^_ exp(-Ef^ (cu(B„)) - Cdm„L) _ 

feasibleexp(-E|'^ {T))\£/fd'n. exp(-C,/m„f) 
exp(-E? (a)(B„))) 

(6.17) = -L-exp(m„(C^f-C^L-C^logK|)). 

Now, if Tmax achieves the maximum of 7t“ (a)(B„)T) over T G then: 

(6.18) ;r^_(®(B„))= Y. 

TGiS'Xn\B„ 

:ft)(B/i)T feasible 
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(6.19) 

( 6 . 20 ) 

( 6 . 21 ) 

( 6 . 22 ) 


< 




= 


< isi/r"- 


^An () "^max ) 

exp(-E^(co(B 

n ) "^max co{dA„))) 

j _ £V| _ 

Hd:da>(dA„) feasible (~®®(*5A„))) 

'^P(-Ef^.(®(Bn)) - QinJ) 
^TGjs'b":[T]^^0^ ^exp(—C,^m„L) 


exp(-Ef (cofB,,))) 

<-^- exp(-m„(Cd£ - CdL - Cdlogj^/j)), 

Bn 


where the square block D-condition has been used in (6.21). Therefore, 


(6.23) 


a-""- < <(a)(B„))Zf„exp(Ef„(a)(B„))) < 


where a := e Qiog|i/|) ^ ^ q ^ (logZg^ — logZg^) —>• 0 (thanks 

to Theorem 3.1), we have obtained (6.13). 

We use (6.13) to represent pressure: 

logZf r logZf 

(6.24) P(4>) = lim = hm / 

|B„| Ii^ooj |B„| 


(6.25) 


= lim 

n^oo 


im [ 

->ooJ 


dv 

0 


-log<,_(at(B„))-E^^(at(B„)) 

|B„| 


dv. 


(Here the second equality comes from the fact that is independent of o, and the 

third from (6.13).) Since v is shift-invariant, it can be checked that: 

-E|(a)(B„)) 


B„| 


(6.26) 

and so we can write: 

(6.27) P(‘l’)= lA^dv- 

It remains to show that: 

/ — IWg 

A„ 


lim [ - 

n^ooj 

P(<I>) = J J 


I 


dv= A^dv, 


lim [ 

n^ocj 


log< (®(B„)) 


|B„| 


dv. 


-log< (cu(B„)) 


BJ 


dv = 


jf^rdv. 


Eix (0 € supp(v) and denote c := Cji{v). We will decompose tt® (a)(B„)) as a product 
of conditional probabilities. By (6.11), for any £ > 0, there exists k := kg so that for 
y,z > Ik, |7ry^j.(tu) — ?f(cu)| < e for all co G supp(v). Eorx G B„_i, we denote B^(x) := 
{yG B„_i :y^ x}. Then, we can decompose tt® (a)(B„)) as: 

(6.29) 7r^_(a)(B„)) = <„(a)(5B„)) J] <, (®W|® (B^WU5B„)) 

(6.30) =<J®(5B„)) Yl ^vix),z{x){(yxico)), 

xGB„_i 


where y(x) := In + x and z(x) := In — x, thanks to the MRP property and stationarity of 
the specification. 












REPRESENTATION AND POLY-TIME APPROXIMATION EOR PRESSURE OE tr LATTICE MODELS 


19 


Let’s denote Rn^k '■= \ ^n-k- Then, B„ = U and we have: 

(6.31) cl^B-.l+l««-Ld Yl < <„(®(B„)) 

(6-32) < Yl 

Taking — log(-), we have that: 

(6.33) 0<-log<(®(B„))- Y. -log7r,w.,w((T.(®)) 

(6.34) <(|5B„| + |7;„_i,,|)log(c-'). 

So, by the choice of k, for a G B„_,(.-b 

(6.35) | 7 t:,(^)_,(^)(( 7 ^(®))-^(c 7 ^(a)))| < e, 

and since 7 t:,.(j.) 2 (j:)(o’j:(®)), 7 f(( 7 i(®)) > c> 0, by the Mean Value Theorem: 

(6.36) |-log7rv(;,),,(;,)(oi(®))-/;r(cJx(a)))| <ec^', 

It follows from (6.11) that n is the uniform limit of continuous functions on supp(v). 
In addition, 7t{co) > c > 0, for all co G supp(v). Therefore, we can integrate with respect 
to V to see that: 


(6.37) 


j - log ((Jx(®))£/v - jl^{co)dv 


< ec 


We now combine the previous equations to see that: 


(6.38) 

(6.39) 


J -\og7:^^{co{B„+i))dv - J Inico)dv\B„_k-i 

< |B„_n|ec-' + (|aB„| + |/?„_i,,|)log (c-i). 


Notice that, for a fixed k, lim„^oo ^ lim„^oo = 1. Therefore, 


(6.40) 

(6.41) 


-ec 


^ + [i 7 i{(o)dv < liminf / 

J n^oo j 

< limsup / 

n—>oo J 

< J Ij:iC0)dV- 


-logKlicoiBn)) 

|B„| 

-log< (®(B„)) 


(6.42) 

By letting e -A 0, we see that: 

-logTT® (®(B„)) 


~ec 


|B«| 

-I 


dv 

dv 


(6.43) lim / 

n^oo J 

completing the proof. 


|B„| 


dv = 


J Ij,{co)dv, 


□ 


7. Spatial mixing and stochastic dominance 

From now on, when talking about specifications for the Potts, Widom-Rowlinson and 
hard-core lattice models, we will distinguish them by the subindex corresponding to the 
parameter j3, A or 7 of the model, i.e. Tip ^ should be understood as a probability measure 

in the Potts model, ^ in the Widom-Rowlinson and in the hard-core lattice model, 
and Tip, Tti and TTy will denote the corresponding specifications. Also, we will write 
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0 and for the functions tta, ft and Ijc in the Potts model, and short-hand notations when 
A = or Sy^^. For example, 

:= Ttlico) := = ®( 0 )) 

The analogous notation will be used for the Widom-Rowlinson and hard-core cases, but 
using the parameters A, and 7 , respectively. 


7.1. Spatial mixing properties. We now introduce concepts of spatial mixing that we 
will need in this paper. Let f : N —>■ R>o be a function such that f(n) \ 0 as n 

Definition 7.1. Given A C we say that a A-MRF ji satisfies strong spatial mixing 
(SSM) with rate /(«) for a class of finite sets if for any A G such that A C A, any 
0 C A, 0 G and ^, t] G with > 0, 

(7.1) |m(0|^)-m( 0|J7)I < |©|/(dist(0,E,A(^,77))). 

We say that a Gibbs specification K = satisfies SSM with rate f{n) for a class 

e 

of finite sets if each element satisfies SSM with rate f{n) for the class ^ restricted to 

subsets of A. 

If there exists C, Ct > 0 such that f can be chosen to be f{n) = Ce^“", we say that 
exponential SSM holds. 


Definition 7.2. ([3, p. 445]) A if-MRF p satisfies the ratio strong mixing property for a 
class of finite sets if there exists C,a > 0 such that for any A G any 0,E C A and 
with p(^)>0, 

piADBl^} 


sup 




- 1 


(7.2) 


: A G .^&,B G .^E,p(Al^)p(Bl^) > 0 
<C ^ 

xe0,yeL 


Proposition 7.1. Let p be a 1?-MRF with supp(/r) = that satisfies the ratio strong 
mixing property for the class of finite simply lattice-connected sets. Then, p satisfies expo¬ 
nential SSM for the family of sets {5y^z}y,z>o. 


Proof. Fix y,z > 0 and the corresponding set Sy^^ g I?. Let 0 C Sy^^, 9 G . 2 /® and ^ 1,^2 € 
with p(^i)p(^ 2 ) > 0 , consider: 

(1) thesetsI:=E 55 ^_^(^i,^ 2 ) and A := 

(2) an arbitrary configuration | G such that | (55j,_^\E) = ^i{dSy^^\'L) (i = 1,2), 
and 

(3) theeventsA := [0] G andB, := [^,(E)] G for i = 1,2. 

Notice that A is a finite simply lattice-connected set and, since supp(/r) = , we can 

be sure that /r(|) >0. Then: 


(7.3) 

(7.4) 

(7.5) 


\pm,)-pm2)\ 


M(A|[|]nBi)-M(A|[|]nB2) 


pjAnBil^) 

M(Bill) 


At(A||)-t-M(A||) 


M(Anfi 2 ll) 


pjADBil^) 


M(AnR 2 ||) 

PiB2\i)pm) 
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(7.6) 

< 2C 


xe&,yeZ 

(7.7) 



rer 


W.l.o.g., we can assume that |E| = 1 (see [9, Corollary 2]). Therefore, by taking C = 2C, 
we have: 

(7.8) |m(0|^i)-ai(0|&)| < 


Remark 2. The proof of Proposition 7.1 seems to require some assumption on the sup¬ 
port of p (for the existence of ^ in the enumerated item list above). Fully supported (i.e. 
supp(/i) = suffices, and is the only case in which we will apply this result (see Corol¬ 
lary 1 ), but the conclusion probably holds under weaker assumptions. 


Given y,z > 0, we define the bottom boundary of ■= ^^y,z © i-®- the 

portion of the boundary of included in the past, and the top boundary of as the 
complement := dSy^^ \ Clearly, dSy^^ = d^Sy^y U 


Proposition 7.2. Let n be a specification satisfying exponential SSM with parameters 
C,a > 0. Then, for all n G N, y,z > In and a G : 


(7.9) 


J/J 


(0(O) = «)-;r,7(0(O) = fl) 


< Ce 


uniformly over 01 , 0)2 € Q.(,(a) such that = 02 (.^). 


Proof. Fix n G N, y,z > In, a G £/ and 01,02 S Q.(,(a) such Oi(.^) = 02 (<^)- Denote 
^ := (Oi{dSn). Then: 


(7.10) 

(7.11) 

(7.12) 

(7.13) 


nl'i0{0)=a)-ngm0)=a) 

= Klie{0)=a)-Y^7t^-^{ei0) = a\ti)7t^^iii) 

< E 4(0(O)=o)-<(0(O)=a) 7r,^^(t7) 

7 ] 


where the summation is taken over all 77 G such that n'^^{r\) > 0 and 77 {d^Sn) = 
O 2 [d^S„y The last inequality abo ve follows from the fact that for any such 77, 5 ^ (^,77) C 

d^Sn, so: 

(7.14) dist( 0 ,E 55 „(^, 77 )) > dist(0,5t5„) =n. 


□ 


Definition 7.3 (Variational distance). Let S be a finite set and letX\ andXz be two S-valued 
random variables with distribution pi and pz, respectively. The variational distance drv of 
Xi andX 2 (or equivalently, of pi and P 2 ) is defined by: 

dTviPuPi) •= 9 E IPiW “P 2 (-)c)| • 

^ xes 


(7.15) 
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It is well-known that dTv{P\ , Pi) is a lower bound on P(Xi ^ X 2 ) over all couplings P 
of Pi and P 2 and that there is a coupling, called the optimal coupling, that achieves this 
lower bound. 

Given a Gibbs specification n, we define: 

( 7 . 16 ) Q{k) := ■ 

The following result is essentially in [7], 

Theorem 7.3. Let 7t be a Gibbs specification for a n.n. interaction <I> and a set of con¬ 
straints S, such that D,{S') has a safe symbol. Then, if pc denotes the critical value of site 
percolation on 1? and Q{Tt) < pc, we have that n satisfies exponential SSM. 

Proof Take p any n.n. Gibbs measure for <I>. Since has a safe symbol, p is fully 

supported, i.e. supp(p) = n(<o’) (very special case of [35, Remark 1.14]). Given a Z'^-MRF 
p, define: 

(7.17) e(p) := maxdrv(p(e(0) = ■lvi),p(9(0) = -1^2)), 

’ii.m 

where rp and p 2 range over all configurations on (9{0} such that p(pi),p(p 2 ) > 0. Then, 
Q(p) < Q(^) < Pc, so by [7, Theorem 1] and shift-invariance of <I>, p satisfies exponential 
SSM (see [30, Theorem 3.10]). Finally, since p is fully supported, we can conclude that n 
satisfies exponential SSM. □ 


7.2. Stochastic dominance. Suppose that ji/ is a finite linearly ordered set. Then for any 
set L (in our context, usually a set of sites or bonds), is equipped with a natural partial 
order which is defined coordinate-wise: for 0i, 02 G we write 0i ^ 02 if 0i(a:) < 
02(x) for every x G L. A function / : —>■ R is said to be increasing if /(0i) < fiOi) 

whenever 0i ^ 82 ■ An event A is said to be increasing if its characteristic function Xa is 
increasing. 

Definition 7.4. Let pi and p 2 be two probability measures on We say that pi is 
stochastically dominated by p 2 , writing pi <0 P 2 , if for every bounded increasing function 
f : > M we have Pi(/) < Piif), where p{f) denotes the expected value Ep(/) of f 

according to the measure p. 


7.2.1. Stochastic dominance and connectivity decay for the bond random-cluster model. 
Recall from Section 5.1 the bond random-cluster model on finite subsets of I? with bound¬ 
ary conditions / = 0,1, and the bond random-cluster model on 1? (see page 11). 


Theorem 7.4 ([17, Equation (29)]). For any p € [0,1] and q GN, and any A C A g 7?: 

(7-18) '^'S.a «««' 41a <d 41a- 

In particular, if p < Pc{q), we have that, for any A d 1?: 

(7-19) 411a<d^p,c,<d411a^ 

where <d A with respect to the restriction of each measure to events on £’**(A). 


The following result was a key element of the proof that fidq) — log(l -f y/q) is the 
critical inverse temperature for the Potts model. We will use this result in a crucial way. 

Recall that for p < pc{q), Pp,ci is the unique bond random cluster measure with param¬ 
eters p and q. 
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Theorem 7.5 ([6, Theorem 2]). Let ^ > 1. For any p < Pc{q) = T+^’ two-point 
connectivity function decays exponentially, i.e. there exist 0 < C{p,q),c{p,q) < °° such 
that for any x,y € 7?: 

(7.20) < C{p,q)e-^^P^‘>'^\\^~y\\f 

where {x O y} is the event that the sites x and y are connected by an open path and || • ||2 
is the Euclidean norm. 


7.2.2. Stochastic dominance for the site random-cluster model. 

Lemma 7.6. Given a set A and parameters p £ [0,1] and q > 0, we have that for 
any X G A and any T £ {0,1}^\W; 

(7.21) pi(q) < ¥p,Ia(^(x) = 1|t) </72(^), 

where pi (q) = and p 2 {q) = pq/ij^py 7« consequence, 

(7.22) <D <D ¥p2iq),A- 
(Recall that 'Lp^A denotes Bernoulli site percolation). 

Proof. This result is obtained by adapting the discussion on [20, p. 339] to the wired site 
random-cluster model. See also [21, Lemma 5.4] for the case q = 2. □ 


7.2.3. Stochastic dominance for the Potts model. As before, \stq£ denote a fixed, but 
arbitrary, choice of a colour. Let A g and consider g : be defined by: 


(7.23) 


{8{mx) = 


if 9(x) =q, 
if 9(x)y^q. 


The function g makes the non-g' colours indistinguishable and gives a reduced model. 
We say 0 ~ 0' if g(9) = g(9'). This relation defines a partition of 32^^ and unions of 
elements of this partition form a sub-algebra of 32^'^, which can be identified with the 
collection of all subsets of {-f, —}^. Let := be the push-forward measure, 

which is nothing more than the restriction (projection) of to {-b, — }^. Chayes showed 
that the FKG property holds on events in this reduced model. In particular: 


Proposition 7.7 ([11, Lemma on p. 211]). For all j5 > 0 and A g I?, satisfies the 
following properties: 

(1) For increasing subsets A, B C {-b,— }^.' 71^ ^(A \ B) > ;r^^(A). 

(2) If A is decreasing and B is increasing, then: 7r^^(A | B) < 7r^^(A). 

(3) If A C A and A is an increasing subset o/{-b, — then: 7r^^(A) > 7r^^(A). 
Proof. 

(1) This is contained in [11, Lemma on p. 211]. 

(2) This is an immediate consequence of (1). 

(3) This is a standard consequence of (1): Let B = -b'^^. Since g^^iji) is a sin¬ 
gle configuration, namely q^^^, we obtain from the Markov property of that 
Tip ^{A) = 7tpj^^{A I B). From (1), we have 7tpj,^(A \ B) > ;r^^(A). Now, combine 
the previous two statements. 


□ 
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Remark 3. The preceding result immediately applies to for events in that are 
measurable with respect to {+, —}^, viewed as a sub-algebra of 

7.2.4. Volume monotonicity for the Widom-Rowlinson model with 2 types. For the classical 
Widom-Rowlinson model {q = 2), Higuchi and Takei showed that the FKG property holds. 
In particular. 

Proposition 7.8 ([21, Lemma 2.3]). Fix q = 2 and let A C A (s and A > 0. Then: 
(7.24) nl^{(Oci)<nl{(Oq). 

However, this kind of stochastic monotonicity can fail for general q (see [17, p. 60]). 

8. Exponential convergence of tt,, in I? lattice models 

In this section, we consider the Potts, Widom-Rowlinson and hard-core models and 
establish exponential convergence results that will lead to pressure representation and ap¬ 
proximation algorithms for these lattice models. 

Recall that for the Potts model, 7r^^(a)) = = (o(0)) and, in particular, 

Ttj^ (co) = = fti(O)), with similar notation for the Widom-Rowlinson and hard 

core models. 


8.1. Exponential convergence in the Potts model. 


Theorem 8.1. For the Potts model with q types and inverse temperature p, there exists a 
critical parameter Pc{q) > 0 such that for 0 < j3 ^ Pc{t]), there exists C,0C > 0 such that, 
for every y,z > In: 


( 8 . 1 ) 




< Ce-“". 


Proof In the supercritical region p > Pc{q), our proof very closely follows [10, Theorem 
3], which treated the Ising case. We fill in some details of their proof, adapting that proof 
in two ways: to a half-plane version of their result (the quantities in (8.1) are effectively 
half-plane quantities) and to the general Potts case. For the subcritical region ^ 
the proposition will follow easily from [3, Theorem 1.8 (ii)]. 


Part I: j3 > Pdq)- Let T^* denote the event that there is a *-path of — from 0 to dSn, 
i.e. a path that runs along ordinary I? bonds and diagonal bonds where the colour at each 
site is not q (in our context below, the configuration on the bottom piece of dS,, will 
be all q and thus a *-path of — from 0 to dS,, cannot terminate on d^Sn). Note that T^* 
is an event that is measurable with respect to the sub-algebra {-b, — }^, for any finite set A 
containing introduced in Section 7.2.3 (recall that this sub-algebra corresponds to the 
reduced Potts model). 

By decomposing TlydfOq) into probabilities conditional on T^^ and (T^* )'^, we obtain: 

( 8 . 2 ) nl{(0ci)-4d(0ci) 

(8.3) = n;^y.^{e{0)=q)-n;^^^^^^{e{0)=q) 

= ("wbofO) = «> - = «ILT.)) 

(8.4) + (1 - (>'?*,('>(«) = «) - = «I(TS)')) . 
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We claim that the expression in (8.2) is nonnegative. To see this, observe that the 
events {w(0) = q\, and may be viewed as the events 

{a)(0) = +}, and {®((95„) = in the sub-algebra of 

the reduced model, as discussed in Section 7.2.3. Now, apply Proposition 7.7 (part 3) and 
Remark 3. 

We next claim that: 


(8.5) 


*,,(0(o)=^i(T.;„r)>*„(0(o)=?). 


To be precise, first observe that co G y iff CO contains an all-^ path in S,, from 
d^n{xi <0} to523^n{2i:i >0}. So, (T^^ )'^ can be decomposed into a disjoint collection 
of events determined by the unique furthest such path from 0. Using the MRF property of 
Gibbs measures, it follows that we can regard each of these events as an increasing event 
in {+, Now, apply Proposition 7.7 and Remark 3. (The reader may notice that here 
we have essentially used the strong Markov property (see [18, p. 1154]).) 

Thus, (8.4) is nonpositive. This, together with the fact that (0(0) = ) = 0, 

yields: 

(8.6) 0 < {(0,) - nUco,) < ;r;^^^^(T-J;r;^^(e(0) = q)< 


•^y,z 


So, it suffices to show that sup^, ht'^p‘'s (^ 5 s ) decays exponentially in n. Fix y,z> In 
and let m > n such that \m >y,z. By Proposition 7.7 (parts 2 and 3) and Remark 3, 


(8.7) 

( 8 . 8 ) 




P:Sy, 


(T^I) < 


= 


13 .B, 






(T, 


dB„) 


So, it suffices to show that sup„,j.„ (T^* ) decays exponentially in n. Recall the 
Edwards-Sokal coupling for the Gibbs distribution and the corresponding bond 

random-cluster measure with wired boundary condition (see Section 5.1). 

W.l.o.g., let’s suppose that n is even, i.e. n =2k < m, for some k gN. We consider the 
following two events in the bond random-cluster model, as in [11, Theorem 3]. Let !k„ be 
the event of an open circuit in B 2 k \ that surrounds B^;.. Let be the event in which 
there is an open path from some site in Bj. to 5B,„. The joint occurrence of these two events 
forces the Potts event (T^* in the coupling: C (T^* )‘^ (here, technically, we 

are identifying these events with their inverse images of the projections in the coupling). 

Then, by the coupling property: 

*,.((T5,*.)')=<‘i,.,.((T5;.)') 

> C.B. ((Ti)'!*" nM„,„) p™ (S. nM.„) 

(8-11) =0i‘i,B„,(3^"nM„,„O, 


so: 


(8-12) ^;k(T,B*„) < i-<(>iIi,B,„(3^«nM„.„) < 

Therefore, 

_(iJcf /r-T'—Jr \ ^ j. ( n 




^(1) 


(8.13) 


(T^b,.) < sup (P^'Ib^ K) + sup 0^;^ (Myj. 
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The first term on the right hand side of (8.13) is bounded from above as follows: 

(S->5) < E 

xGdBi^,y€:^B2i( 

(8-16) < Y, ^P’’,qix^y)y 

x€dBii,y€^B2k 

where Bm = [—m+ l,m]^ and p* denotes the dual of p and the inequalities follow 
from Proposition 5.3 and Theorem 7.4. 

If p > Pc{q), then p* < Pc{q), and by Theorem 7.5, the first term on the right side of 
(8.13) is upper bounded by 64C(p*,g')«^exp(—c(/7*,^)«/4), since |<3 B;i|| 5B2A:| < 64n^ and 
\\x — y\\ 2 >k—\ > for all x G ^Bj. and y € ^B 2 /c. So, the first term on the right side of 
(8.13) decays exponentially. 

As for the second term, in order for M„ to fail to occur, there must be a closed circuit 
in Bm \ B|(. and in particular a closed path from L,„.„ := B^ \ B,;. fl {xi < 0 ,X 2 = 0} to 
Rm,n := B„, \ B^. n {xi > 0,X2 = 0} in Bm- Thus, 


(8.17) 


(8.18) 

w 

VI 


xeL,„^„,yeR, 

(8.19) 

w 

VI 




,9,B„,+1 ' 


where the last inequality follows by Proposition 5.3 and Proposition 7.4. By Theorem 7.5, 
this is less than: 


( 8 . 20 ) 

( 8 . 21 ) 


Y C{p* 


_ ^{P —2c(p*,q)n 


2 


Thus, the 2nd term on the right side of (8.13) decays exponentially, so sup„j>„ ) 

decays exponentially in n. Thus, by (8.7) sup,„>„ decays exponentially in 

n, as desired. 


Part II: j3 < Pdq)- Recall from Section 7 the notions of strong spatial mixing and ratio 
strong mixing property. 

Theorem 8.2 ([3, Theorem 1.8 (ii)]). For the 1? Potts model with q types and inverse tem¬ 
perature p, ifO < P < Pc{q) o.nd exponential decay of the two-point connectivity function 
holds for the corresponding random-cluster model, then the (unique) Potts Gibbs measure 
satisfies the ratio strong mixing property for the class of finite simply lattice-connected 
sets. 

Corollary I. For the 1? Potts model with q types and inverse temperature 0 < P < Pc{q), 
the specification satisfies exponential SSMfor the family of sets 

Proof This follows immediately from Theorem 7.5, Theorem 8.2 and Proposition 7.1. □ 
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Then, since exponential SSM holds for the class of finite simply lattice-connected sets 
when p < Pc{q), the desired result follows directly from Proposition 7.2. 

This completes the proof of Theorem 8.1. □ 


8.2. Exponential convergence in the Widom-Rowlinson model. Recall that for Bernoulli 
site percolation in 1? there exists a probability parameter pc, known as the percolation 
threshold, such that for p < pc, there is no infinite cluster of Ts 22 ‘^ltnost surely and 
for p > Pc, there is such a cluster xj/p ^z-almost surely. Similarly, one can define an analo¬ 
gous parameter p* for the lattice 7?’*, which satisfies Pc + Pc = ^ (see [37]). 


Theorem 8.3. For the Widom-Rowlinson model with q types and activity X, there exist two 
critical parameters 0 < Ai {q) < X 2 {q) such that for A < Ai {q) or X > X 2 {q), there exists 
C,a > 0 such that, for every y,z> In: 


( 8 . 22 ) 




< Ce-“". 


Proof As in the proof of Theorem 8.1, we split the proof in two parts. 


Part I: A > X 2 {q) := q^ j • « € N andy,z > In. Notice that, due to the constraints 

of the Widom-Rowlinson model, and recalling Proposition 5.4: 

(8.23) =q) = = 1), 

where p = and the same holds for n^{o}q). Then, it suffices to prove that: 


(8.24) 


v^J!k(0(O) = i)-v^J!k,(0(O) = i) 


< Ce 


for some C, a > 0. 

Notice that 0 G S„ ^ =: A. Fix any ordering on the set A. From now on, when we 

talk about comparing sites in A, it is assumed we are speaking of this ordering. For con¬ 
venience, we will extend configurations on S„ and A to configurations on A by appending 
lA\ 5 „ respectively. 

Now, we will proceed to define a coupling lPn,y,z of j and defined on pairs 

of configurations (0i, 62 ) G {0,1}^ x {0, l}''^. The coupling is defined one site at a time, 
using values from previously defined sites. 

We use (t{ , T 2 ) to denote the (incomplete) configurations on A x A at step f = 0,1,..., |. 

We therefore begin with tJ* = 1'^^'^" and T 2 = Next, we set Tj = rj* and form Z 2 by 
extending T 2 to A \ Sn, choosing randomly according to the distribution A (’ I ^ 

this point of the construction, both Tj and rj have shape A\5„. In the end, 
will give as a result a pair ( 0 i, 02 )- 

At any step t, we use W' to denote the set of sites in A on which z\ and Zj have already 
received values in previous steps. In particular, W* = A\5„. At an arbitrary step t of the 
construction, we choose the next site on which to assign values in and as 
follows: 

(i) If possible, take to be the smallest site in d*W^ that is ^-adjacent to a site 
y GW’ for which {T:\iy), zl^iy)) ^( 1 , 1 ). 

(ii) Otherwise, just take to be the smallest site in d*W’. 

Notice that at any step t, W‘ is a ^-connected set, and that it it always possible to find 
the next site for any t < |5„| (i.e. the two rules above give a well defined procedure). 
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Now we are ready to augment the coupling from to U by assigning 




^pisyX- I ^ 2 ) 


*) according to an optimal coupling of (’KS) 


{^+ 1 } 


and 


) and T 2 *(-r‘ 

, i.e. a coupling which minimizes the probability that, given (t{,T 2 ), 
'n,y^^ is defined site-wise, and at each step is assigned accord- 


0 i(x*+i) ^ 02 (x'+*). Since 


( 1 ) 


x\) in the first coordinate and ^ (• | T 0 in the second, the reader may 


ing to i-rv .- --- — 

check that it is indeed a coupling of 1 //^*^ and ^ 
following. 


. The key property of 1 


'hy,z 


is the 


Lemma 8.4. 0i(O) ^ 02 ( 0 ) if and only if there exists a path T of -k-adjacent 

sites from 0 to dSn, such that for each site y €T, ( 0 i (y), 02 (y)) f ( 1 , 1 )- 


Proof Suppose, for a contradiction, that 0i (0) f 02 (0) and there exists no such path. This 
implies that there exists a circuit C surrounding 0 (when we include the bottom boundary 
as part of C) and contained in Sn such that for all y £ C, (0i (y), 02 (y)) = (1 1 1 )■ Define by 
I the simply lattice-*-connected set of sites in the interior of C and, let’s say that at time to, 
x'o was the first site within I defined according to the site-by-site evolution of T‘n,y,z- Then, 
( (y®), t 2 )) cannot have been defined according to rule (i) since all sites *-adj acent to 

x'® are either in I (and therefore not yet defined by definition of x^®), or on C (and therefore 
either not yet defined or sites at which 0 i and 02 are both 1 ). 

Therefore, (0i (x'®), 02 (x:'®)) was defined according to rule (ii). We therefore define the 
setD;=A\W®^' D /, and note that 0 and x*® belong to the same ^-connected component 
0 of D. We also know that t[“ * (<3*D) = 1'^*^, otherwise some unassigned 

site in D would be *-adjacent to a 0 in either rj” *((9*D) or x'f ^{d*D), and so rule (i) 
would be applied instead. We may now apply Lemma 5.5 (combined with Remark 1) to 
© and A in order to see that (0i(®)Ki” *) ( 02 (®)|'f 2 ” *) identical. 

This means that the optimal coupling according to which xf (x'®) and X 2 (x^®) are assigned 
is supported on the diagonal, and so t|'’(x:'®) = T 2 ’(x:'®), Pn^y^j-almost surely. This will 
not change the conditions under which we applied Lemma 5.5, and so inductively, the 
same will be true for each site in I as it is assigned, including 0. We have shown that 
01 ( 0 ) = 02 ( 0 ), Pn.j, ,-almost surely, regardless of when 0 is assigned in the site-by-site 
evolution of Vn,y,z- This is a contradiction, and so our original assumption was incorrect, 
implying that the desired path T exists. □ 


Given an arbitrary time f, let: 


(8.25) 


pS(-):= V/'S,5„(- 




{.T^} 


) r..n Pli-) ■= 


{■v'j 


be the two corresponding probability measures defined on the set {0,1}^’^^. Note that 
at any step within the site-by-site definition of Pn^,;, Lemma 7.6 implies that < 
p)(l), where X = and i = 1,2. Now, w.l.o.g., suppose that P 2 ( 0 ) > p5(0). Then, 
an optimal coupling Q' of pj and Pj will assign Q'({( 0 , 0 )}) = p 5 ( 0 ), Q'({( 0 , 1 )}) = 0 , 
Q'({(li0)}) =P 2 ( 0 )-Pi( 0 )’ and Q'({(1)1)}) = 1 -P 2 (oi Therefore, 


X+q^' 


(8.26) 


Q'({(l,l)r)=P2(0)< 








REPRESENTATION AND POLY-TIME APPROXIMATION EOR PRESSURE OE tr LATTICE MODELS 


29 


Next, define the map h ; {0,1}'^" x {0,1}^" —?► {0, l}*^" given by: 

(8.27, (/,«)„9.))(.v) = |‘ 

|0 if((,,(i),(,2(»)) = (l,l). 


By (8.26), (the push-forward measure) can be coupled against an i.i.d. measure 

on { 0 , 1 }^" which assigns 1 with probability and 0 with probability and that 

the former is stochastically dominated by the latter. This, together with Lemma 8.4, yields 


(8.28) 

(8.29) 


= 1) - = 1)| < P«.rx(0i(O) 02(0)) 


Since we have assumed X> Pc +/ 2 j = B we have < p* p follows 

by [1, 32] that the expression in (8.29) decays exponentially in n. This completes the proof. 


Part II: X < Xi (q) := ^ j • Observe that, by virtue of Proposition 7.2, it suffices to 

prove that tt™ satisfies exponential SSM. For this, we use Theorem 7.3. By considering 
all cases of nearest-neighbour configurations at the origin, one can compute: 


(8.30) 


e(^r) 


max dTv{^^' 
mi,o>2eQ.{S‘) 


m \ 
'A.{0}’"2.,{0}^ 


qX 

1 -|- qX 


By Theorem 7.3, we obtain exponential SSM when: 

(8.31) x<-(-^\=X,{q). 

q\^-pcj 

Uniqueness of Gibbs states in this same region was mentioned in [19, p. 40], by appeal¬ 
ing to [7, Theorem 1] (which is the crux of Theorem 7.3). □ 


Remark 4. In the case q = 2, it is possible to give an alternative proof of Theorem 8.3, 
Part I, using the framework of the proof of Theorem 8.1, Part I. The arguments through 
(8.7) go through, with an appropriate re-definition of events and use of Proposition 7.8 
for stochastic dominance. One can then apply Lemma 7.6 to give estimates based on the 
site random-cluster model. (In contrast to Theorem 8.1, Part I, this does not require the 
use of planar duality). So far, this approach is limited to q = 2 because we do not know 
appropriate versions of Proposition 7.8 for q > 2. 


8.3. Exponential convergence in the hard-core model. Our argument again relies on 
proving exponential convergence for conditional measures with respect to certain “ex¬ 
tremal” boundaries on S„, but these now will consist of alternating 0 and 1 symbols rather 
than a single symbol (recall from Section 4.3 that 0 )^"^ is defined as the configuration of 
I’s on all even sites and 0 on all odd sites). 


Theorem 8.5. For the 1? hard-core model with activity y, there exist two critical param¬ 
eters 0 < 7 i < 72 such that for any 0 < 7 < 71 ory>yi, there exist C, a > 0 such that for 
every y,z> In, 


< Ce 


— an 


(8.32) 
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Proof. As in the previous two theorems, we consider two cases. 


Part I: 7 > 72 := 468. Our proof essentially combines the disagreement percolation tech¬ 
niques of [7] and the proof of non-uniqueness of equilibrium state for the hard-core model 
due to Dobrushin (see [13]). We need enough details not technically contained in either 
proof that we present a mostly self-contained argument here. From [7, Theorem 1] and 
an averaging argument (as in the proof of Proposition 7.2) on d^Sn induced by a boundary 
condition on Sy_^, we know that for any y,z> In, 


(8.33) 




< ^n,y,z (3 a path of disagreement from 0 ^ 


\y,z of and 


. We do not need the structure of 1 


s„ 


n,y,z 


for a certain coupling 

here, but instead note the following: a path of disagreement for the boundaries co^°\dSn) 
and (dSy^z) implies that in one of the configurations, all entries on the path will be “out 
of phase” with respect to C0^°\ i.e. that all entries along the path will have 1 at every odd 
site and 0 at every even site rather than the opposite alternating pattern of Then, if 
we denote by the event that there is a path T from 0 •(-> d^S„ with 1 at every odd site 
and 0 at every even site, it is clear that: 


(8.34) ^n,y,z (3 a path of disagreement from 0 UA d^S„) < {^n) + ^y^SyA^n)- 

Since y,z> In are arbitrary (in particular, y and z can be chosen to be In), it suffices to 
prove that sup,.->j„ decays exponentially with n. Define the set: 


(8.35) 0v,z = {0 G {0,1}^- : 0 is feasible and e{d*Sy^^) = co^‘’\d*Sy^,)}. 

For any 0 G 0y,z, we define Eo(0) to be the connected component of Ejj,^(0,0)^“^) 
(= {x G : 0(x) f ta^")(x)}) containing the origin 0. Since ^ C {Eq(0) H 0}, 

our proof will then be complete if we can show that there exist C, a > 0 so that for any n 
and 7 ,z > In, the following holds: 


(8.36) (Eo(0) n dySn f%)< Cf-"". 

To prove this, we use a Peierls argument, similar to [13]. 

Fix any 7 ,z > In and for any 0 G 0y,z, define Eo(0) as above, and let fG(0) to be the 
connected component of {x G Sy^ : 9(x) = containing d*Sy^z- Clearly, Eo(0) 

and K{6) are disjoint, K{9) f 0 and, provided 0(0) = 0, Eo(0) f 0. Then, define r(0) := 
T,o{9)r\dK{9) C Sy -. We note that for any 0 G 0v- with 0(0) = 0, we have that 0(r(0)) = 
or(0)^ 

since adjacent sites in Eo(0) and K{9) must have the same letter by definition of 
Eo(0), and adjacent 1 symbols are forbidden in the hard-core model. Therefore, every 
X G r(0) is even. 

We need the concept of inner external boundary for a connected set E (s Z^. The inner 
external boundary of E is defined to be the inner boundary of the simply lattice-connected 
set consisting of the union of E and the union of all the finite components of I? \ E. Intu¬ 
itively, the inner external boundary of E is the inner boundary of the set E obtained after 
“filling in the holes” of E. Notice that the set r(0) corresponds exactly to the inner external 
boundary of Eo(0). In addition, by [12, Lemma 2.1 (i)], we know that the inner external 
boundary of a finite connected set (more generally a finite *-connected set) is *-connected. 
Thus, r(0) C Sy^z is a *-connected set C* that consists only of even sites and contains the 
origin 0, for any 0 G 0y,z with 0 (0) = 0. 





REPRESENTATION AND POLY-TIME APPROXIMATION EOR PRESSURE OE tr LATTICE MODELS 


31 


Then, for C* C 5^, we define the event Eq* := {0 g 0^^ : r(0) = C*}, and will bound 
from above for every C* such that Eq* is nonempty. We make some more 

notation; for every such a set C*, define 0(C*) (for ‘outside’) as the connected component 
of (C*)‘^ containing d*Sy^-, and define /(C*) (for ‘inside’) as \ (C* U (9(C*)). Then C*, 
/(C*), and 0{C*) form a partition of Sy-. We note that there cannot be a pair of adjacent 
sites from /(C*) and (9(C*) respectively, since they would then be in the same connected 
component of (C*)^. We also note that for every 0 g Eq*, C* C Eo( 0) ^ C*U/(C*) 
and K{6) C 0(C*) though the sets need not be equal, since Eo(0) or K(9) could contain 
“holes” which are “filled in” in /(C*) and 0{C*), respectively. 



Figure 2. A configuration 0 g Eq*. On the left, the associated sets 
Eo(0) and K{e). On the right, the sets /(C*) and 0(C*) for r(0) = C*. 


Choose any set C* such that Ec* ^ 0. For each 0 g Eq* and x g C*, using the definition 
ofC* and the fact that/r( 0) C 0{C*), there exists xq g {ei,—ei,e 2 ,—^ 2 } forwhichx —xq g 
(9(C*). Fix an xq which is associated to at least |C*|/4 of the sites in C* in this way. Then, 

we define a function s : Eq* {0, 1 }'^az that, given 0 g Eq*, defines a new configuration 
s(0) as follows; 


(8.37) 


(i(0))(x) 


{ 0(x — Xq) ifxg/(C*), 

0(x) ifxge)(C*), 

1 ifxg C* andx —Xq g (9(C*), 

0 if X g C* and X — Xq g/(C*). 


Informally, we move all 1 symbols inside /(C*) in the XQ-direction by 1 unit (even if 
those symbols were not part of Eo(0)), add new 1 symbols at some sites in C*, and leave 
everything in 0{C*) unchanged. 

It should be clear that s(0) has at least |C*|/4 more 1 symbols than 0 did. We make the 
following two claims; s is injective on Eq*, and for every 0 g Eq*, s{9) g If these 
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claims are true, then clearly (^(iic*)) > 7^ (^c*), implying that: 


(8.38) 



Firstly, we show that s is injective. Suppose that 9i ^ 92, for 0i, 02 G Eq*- Then there is 
a site X at which 0i(x) ^ 02(x). Ifx G 0{C*), then (i(0i))(x) = 0i(x) ^ 02(x) = (i(02))(x) 
and soi(0i) ^i(02). IfxG/(C*), then (i(0i))(x + xo) = 0i(x) 02(x) = (i(02))(x+xo), 

and again s(0i) ^ i(02). Finally, we note that x cannot be in C*, since at all sites in C*, 
both 01 and 02 must have 0 symbols. 

Secondly, we show that for any 0 G E( 2 *, s{9) is feasible. All that must be shown is 
that s(0) does not contain adjacent 1 symbols. We break 1 symbols in s{9) into three 
categories: 

(1) shifted, meaning that the 1 symbol came from shifting a 1 symbol at a site in /(C*) 
in the XQ-direction, 

(2) new, meaning that the 1 symbol was placed at a site x G C* such that x — xq G 
e)(C*), or 

(3) untouched, meaning that the 1 symbol was at a site in 0{C*) (3 d*Sy^f). 

Note that untouched 1 symbols cannot be adjacent to C*: 0 contains all 0 symbols on 
C*, and so since C* C Eo(0), a 1 symbol adjacent to a symbol in C* would be in Eo(0) as 
well, a contradiction since T.o{9) C C* U/(C*), and so Eo(0) and 0(C*) are disjoint. 

Clearly shifted 1 symbols cannot be adjacent to each other, since there were no adjacent 
1 symbols in 0. All new I’s were placed at sites in C*, and all sites in C* are even, so new 
1 symbols can’t be adjacent to each other. Untouched I’s can’t be adjacent for the same 
reason as shifted I’s. We now address the possibility of adjacent 1 symbols in s(9) from 
different categories. A shifted or new 1 in s(9) is at a site in C* U/(C*), and an untouched 
1 can’t be adjacent to a site in C* as explained above, and also cannot be adjacent to a site 
in /(C*) since /(C*) and (9(C’^) do not contain adjacent sites. Therefore, shifted or new I’s 
can’t be adjacent to untouched I’s. The only remaining case which we need to rule out is a 
new 1 adjacent to a shifted 1. Suppose that (s(0))(x) is a new 1 and (s(0))(x') is a shifted 
1. Then by definition, x' — xq G /(C*) andx — xq G 0(C*). We know that /(C*) and 0(C*) 
do not contain adjacent sites, so x — xq and x' — xq are not adjacent, implying that x and x' 
are not adjacent. We’ve then shown that s(9) is feasible and then, since d*Sy^- C 0(C*), 
s(9) G 0y,;, completing the proof of (8.38). 

Recall that every set C* which we are considering is ^-connected, occupies only even 
sites, and contains the origin 0. Then, given k G N, it is direct to see that the number of 
such C* with |C*| = k is less than or equal to k-t{k), where t{k) denotes the number of 
site animals (see [25] for the definition) of size k (the first k factor comes from the fact that 
site animals are defined up to translation, and here given a site animal of size k, exactly 
k translations of it will contain the origin 0). We know that for every e > 0 there exists 
Cg > 0 such that t{k) < CgjS + e)*’ for every k, where d := lim/.^oo (f< 4.649551 
(see [25]). 

If Eo(®) fl ^o(0) to intersect the left, top or right boundary of 5„. 

W.l.o.g., we may assume that Eo(0) intersects the right boundary of 5„. Then, every verti¬ 
cal segment in the right half of must intersect Eo(0) therefore, at least one element 
of its inner external boundary, namely r(0). Then: 


(8.39) 


Io(0)n5|5„ ^0 |r(0)|>n. 
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Therefore, taking an arbitrary e > 0, we may bound (Eo(0) ^ 0) from 

above; 

oo 

(8.40) 7r“J’^(Eo(0)nat5„?^0)< L + £)'•• 7-"^^ 

C*:|C*|>n k=n 

which decays exponentially in n as long as 7 > (5 -f e)^, independently of 7 and z. Since e 
was arbitrary, 7 > 468 > 5^ suffices for justifying (8.36), completing the proof. 

Part II: 7 < 71 := 2.48. It is known (see [39]) that when d = 2 and 7 < 2.48, satisfies 

exponential SSM. Then, by applying Proposition 7.2, we conclude. □ 

9. POLY-TIME APPROXIMATION FOR PRESSURE OF LATTICE MODELS 

By a poly-time approximation algorithm to compute a number r, we mean an algorithm 
that, given S N, produces an estimate rN such that \r — ri^\ <h and the time to compute 
is polynomial in N. 

Theorem 9.1. Let ^ be a n.n. interaction for a set of restrictions S' and suppose that 
Ll{S) satisfies the square block D-condition. Let W € Ll{S) be a periodic point such that 
C 7 r(v®) > 0. In addition, suppose that there exists C,(X > 0 such that, for every y,z > In: 

(9.1) \n„{(o) - ny,fi(o)\ <Ce^“" over co eO{W). 

Then: 

(9.2) P(4>) = i^ E 4M+A<i,(ai), 

and there is a poly-time approximation algorithm to compute P(‘P), when d = 2. 

Proof Notice that supp(v®) = 0(®) (LLl{S), since Q.{S) is shift-invariant and m S Q.{S). 
Now, since |7r„(a)) — 7rv,^(fii)| < Ce^“" over O) € supp(v“), we can easily conclude that 
limy^^^oo ny-{(o) = 7t{co) uniformly over co S supp(v“). This, combined with Q(S) satis¬ 
fying the square block D-condition and > 0, gives us 

(9.3, P(4.)=/(4+^,).v-=p^ I 

' C0€supp(v^) 

thanks to Theorem 6.3. 

For the algorithm, it suffices to show that there is a poly-time algorithm to compute 
^( 01 ), for any co € 0 (m). 

By Equation 9.1, there exist C,a > 0 such that \7i„{co) — 7 r(fti)| < Since |(95„| 

is linear in n when if = 2, by a modified transfer matrix approach (see [31, Lemma 4.8]), 
we can compute 7tn{co) in exponential time Ke^" for some K,p >0. Combining the ex¬ 
ponential time to compute 7tn{co) for the exponential decay of |7r„(a)) — 7t{co)\, we get a 
poly-time algorithm to compute P(<I>): namely, given N € N, let n be the smallest integer 
such that Then 7r„+i(a)) is within ^ of Tt{co) and since ^ < Ce^“", the 

time to compute 7r„+i (w) is at most; 

(9.4) = {KePCPl^)- — ^ ^ , < {KeP!^)NP>^, 

(Ce 


which is a polynomial in N. 
Corollary 2. The following holds: 


□ 
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(1) For the I? Potts model with q types and inverse temperature j3 > 0 .- 

(9.5) P(<J>^)=/f(®,)+2j3. 

(2) For the 1? Widom-Rowlinson model with q types and activity A £ (0,Ai(q'))U 

(9.6) P(O0=/x(®./)+log^: 

where Xi{q) := | (t^) and Xjiq) ■= q^ (t^)' 

(3) For the 1? hard-core model with activity y £ (0, 7 i) U ( 72 ,°°).' 

(9.7) p(<i>^)^l/^(a)W) + llog 7 , 
where 71 = 2.48 and 72 = 468. 

Moreover, for the three models in the corresponding regions (except in the case when 
j3 = j3e(q') in the Potts model), the pressure can be approximated in poly-time, where the 
polynomial involved depends on the parameters of the models. 

Proof The representation of the pressure given in the previous statement for the 7? Potts 
model with q types and inverse temperature )3 7 ^ ^c{q), the I? Widom-Rowlinson model 
with q types and activity X £ (0, Ai(q'))U (A 2 (q'),°°) and the 1? hard-core model with activ¬ 
ity 7 G (0, 7 i) U ( 72 ,°°), is a direct consequence of Theorem 9.1, by virtue of the following 
facts: 

• Recall that the corresponding n.n. SFT for the Potts, Widom-Rowlinson 

and hard-core model has a safe symbol, respectively, so il(tf) satisfies the square 
block D-condition and Cn:(v) > 0, for any shift-invariant v with supp(v) C 

in each case. 

• If we consider the delta-measure v = v“'> = 5(b,, both in the Potts and Widom- 

Rowlinson cases (in a slight abuse of notation, since the Potts and Widom-Row¬ 
linson a-algebras are defined in different alphabets), or the measure v = = 

^S^(e) -h yS^(o) in the hard-core case, we have that in all three models, for the 
range of parameters specified, except for when j3 = /3c (q) in the Potts model, there 
exists C,a > 0 such that, for every y,z> in: 

(9.8) | 7 r„(a)) — 7 ry,j(tu)| < over to G supp(v), 

thanks to Theorem 8.1, Theorem 8.3 and Theorem 8.5, respectively. (Notice that 
/J(®W)=A4.(®W)=0.) 

This proves (9.5), (9.6) and (9.7), except in the Potts case when j3 = j3c. To establish 
this case, first note that it is easy to prove that P(<t>j 3 ) is continuous with respect to j3. 

Second, if j3i < j 32 , then (o^) < This follows by the Edwards-Sokal coupling 

(see Theorem 5.2) and the comparison inequalities for the bond random-cluster model [2, 
Theorem 4.1]. 

As an exercise in analysis, it is not difficult to prove that if am,n > 0, and each am+i,n < 
affi.fi and a^ji a^yi^yi, then liiU/ji lim^j a/^i^fj — lintj] lim^ ayyi^fi — a, for some a 0 . 

Now, consider the sequence := (cu^). By stochastic dominance (see Propo¬ 

sition 7.7), am,n is decreasing in n. By the previous discussion (Edwards-Sokal coupling), 
it is also decreasing in m. Therefore, and since am,n > 0, we conclude that lim,„ lim„ a,n.n = 
lim„ limrn am,n = a, for some a. 
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Then, we have that; 

(9.10) =lim-loglim;r,f"*^^^^”(a)<y)+2 (Pciq) + — 

m n \ m 


(9.11) 

= -loglimlim;r,f"*''^^^"' {cOq)+2^c{q) 

m n 

(9.12) 

= -loglimlimTT^"*'^^^"' {(Oq)-\-2j5c{q) 

n m 

(9.13) 

= -\og\m\nt^‘‘\(Oq) -y2^^{q) 

n 

(9.14) 

= lt^‘^\c0q)+2^(q). 


(To prove that linim {(Oq) = ni’'^‘‘\(Oq) is straightforward.) 

Finally, the algorithmic implications are also a direct application of Theorem 9.1. □ 

Remark 5. The algorithm given in Theorem 9.1 seems to require explicit bounds on the 
constants C and a, so that given N € N, we can find an explicit n such that < jj. 

Without such bounds, while there exists a poly-time algorithm, we do not always know how 
to exhibit an explicit algorithm. However, for all three models, for regions sufficiently deep 
within the supercritical region (i.e. j3, X or y sufficiently large), one can find crude, but 
adequate, estimates on C and a and thus can exhibit a poly-time algorithm. This is the case 
for the hard-core model, where our proof does allow an explicit estimate of the constants 
for any y > 468. On the other hand, in the regions specified in Corollary 2 within the 
subcritical region, all three models satisfy exponential SSM and then using [31, Corollary 
4.7], one can, in principle, exhibit a poly-time algorithm (even without estimates on C and 
a). 
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